Frustrated kinetic energy, the optical sum rule, and the mechanism of 

superconductivity 
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The theory that the change of the electronic kinetic energy in a direction perpendicular to the CuO- 
planes in high-temperature superconductors is a substantial fraction of the condensation energy is 
examined. It is argued that the consequences of this theory based on a rigorous c-axis conductivity 
sum rule are consistent with recent optical and penetration depth measurements. 



PACS: 

a. Introduction. The aim of this Letter is to partly 
resolve a number of issues related to a theory of 

high-temperature superconductivity known as the inter- 
layer tunneling theory (ILT) Jl(J and to propose the ef- 
ficacy of a conductivity sum rule. Within a simple ver- 
sion of ILT, one relates the zero-temperature c-axis pen- 
etration depth A c to the superconducting condensation 
energy H. Here, we point out that the realization of 
ILT and the interpretation of recent measurements of A c 
|^|,^), necessarily require more careful analysis and that 
the two can be brought into agreement. In addition, we 
argue that ILT accounts for two features of c-axis op- 
tical measurements: (1) the observation that the c-axis 
(perpendicular to the CuO-planes) kinetic energy is sub- 
stantially reduced in the superconducting state B and 
(2) the correlation ("Basov correlation") between A c and 
the c-axis conductivity in the normal state pd[ | . 

For our purposes, the content of ILT is that a signifi- 
cant portion of the superconducting condensation energy 
comes from the change in the c-axis kinetic energy as the 
electrons enter the superconducting state. It is a phe- 
nomenological fact that this kinetic energy is frustrated 
in the normal state, but that the frustration is relieved 
in the superconducting state as the coherent tunneling of 
pairs becomes possible, resulting in a sharp plasma edge 
in the reflectivity E] , 

Recently, Anderson |J conjectured that the full con- 
densation energy is derived from the c-axis Josephson en- 
ergy, which, in turn, determines the penetration depth. 
Then, using estimates of the condensation energies, he 
predicted A c . On the basis of recent experiments it 
has been suggested that this prediction is strongly vio- 
lated in both Tl 2201 and Hg 1201, although it appears 
to hold for LSCO for a large range of doping. The single- 
layer superconductors containing one CuO-plane per unit 
cell are emphasized because they pose the most stringent 
test of ILT. 

However, the situation is not so clear. (1) As shown 
earlier H, the predicted A c should be a factor 2 larger 
than that predicted in H. (2) The measured values of 
A c are in disagreement. Vortex imaging measurements 



H give A c = 17 - 19 /j,m in Tl 2201, while it is 12 fim 
in the optical measurements in a similar sample. For 
Hg 1201, vortex measurements give 8 /jm ||, while the 
optical measurements give 6.2 /im, again in a similar sam- 
ple E3I, and, disturbingly, magnetic measurements yield 
1.36 fim JlJ]. (3) The normal state electronic specific 
heat must be extrapolated to T = from above T c to 
determine the condensation energy. 

There is an even more fundamental difficulty. The con- 
densation energy is well-defined only within mean field 
theory. For those materials that deviate from mean field 
behavior, that is, those that do not have a sharp spe- 
cific heat jump at T c , the condensation energy cannot be 
determined by a simple integration of the specific heat 
0. The electronic specific heat in Tl 2201 [§J| and in 
Hg 1201 [^6| show large fluctuation effects; we shall show 
that in these cases, agreement with experiments can be 
achieved if the fluctuation effects are subtracted out. 

b. Sum rule. We now discuss the c-axis conductivity 
sum rule J7J. Consider the full hamiltonian H = H rcst + 
H c ; the c-axis kinetic energy is defined by 



He = - 



]t±(ij, l)ct s c ji+hs +h.c. 



(1) 



the remainder, H rcstl contains no interplane interaction 
terms |lj , but it is otherwise arbitrary and may contain 
impurity interactions that couple to the charge density. 
The hopping matrix element t±(ij,l), where refer 
to the sites of the two-dimensional lattice, and I to the 
layer index, can be random in the presence of impurities 
fl8f . The electron operators (c, c^ ) are also labeled by a 
spin index s. We denote the magnitude of t±(ij,l) by 
t±. It is easy to adapt a sum rule derived first by Kubo 
|]l9f to get a sum rule for the c-axis optical conductivity 
a c (uj, T), which is 



du Re cr c {u,T) 



2h z {Ad) 



'-He)- 



(2) 



Here the average refers to the quantum statistical aver- 
age, A to the two-dimensional area, and d to the separa- 
tion between the CuO planes. 
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The hamiltonian H c is an effective hamiltonian valid 
for low energy processes that do not involve interband 
transitions. It can be derived by a downfolding process, 
in which all the higher energy bands are integrated out 
pjj . Because interband processes involve large energy 
differences, a second order downfolding procedure is suf- 
ficient. This is essentially how one derives the effective 
mass contribution of higher energy bands. 

Since H c is a low-energy effective hamiltonian, the up- 
per limit in Eq. (0) can not exceed an interband cut- 
off ui c , of order 2-3 eV. In the superconducting state, 
a cs (uj,T) = D c {T)S(to) + of* g (w, T), where D C (T) is the 
superfluid weight. From Eq. (g), it follows that 



£>c(2i 



+ 



o+ 
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Rea cn (^,T 2 )-Re< c s g (w,T 1 ) 
{-H C {T X )) S - (-H C (T 2 )), 



(3) 



Here, if T 2 < T c , {—H c (T2)) n is to be understood as taken 
in the normal state extrapolated to below T c , and a cn is 
the corresponding conductivity. 

c. Lowering of the kinetic energy in the superconduct- 
ing state: Basov et al. § have tested the sum rule, 
Eq. (|), by setting T x < T c and T 2 = T c . The result 
is that up to 0.15 eV, the integral over the conductivity 
is only half the left hand side, so that the remaining half 
must then be the change in the c-axis kinetic energy of 
electrons. Since we expect the normal state kinetic en- 
ergy to become less negative as the temperature is low- 
ered pc| , this measured change of the kinetic energy is 
only a lower bound on the difference in the kinetic energy 
at T — 0. Thus, these experiments support the funda- 
mental statement of ILT that the c-axis kinetic energy 
is substantially lowered in the superconducting state, in 
contrast to BCS. 

d. The T — superfluid weight D c and X c : We relate 
A c to the change in c-axis kinetic energy as follows: We 
set T\ = T 2 = in Eq. (||). From the experiments of 
Ando et al. |2(J , it is seen that the c-axis response in the 
normal state obtained by destroying superconductivity is 
insulating as T -> 0; it follows that a cn (uj,T = 0) - w 2 , 
as uj — > 0. The regular part a^ g (cu,T = Q) is also ex- 
pected to vanish as a power law in a d-wave supercon- 
ductor P7j|. At high frequencies the two conductivities 
must, however, approach each other. Consequently it 
is reasonable to hypothesize that the conductivity inte- 
gral on the right hand side of Eq. (||) is negligibly small. 
Therefore, 
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where we assumed local London electrodynamics. We 
emphasize that the choice of the normal state in Eq. (0) is 
not arbitrary because we have assumed that the integral 



on the right hand side of Eq. (||) is vanishingly small, 
and this would not be true for an arbitrary state. In any 
case, the right hand side should be a lower bound. 

e. Condensation energy: The attempt to extract the 
condensation energy from the specific heat data runs into 
ambiguity, except within a mean field treatment. In the 
presence of fluctuations, superconducting correlations, 
which can primarily be of in-plane origin, contribute to 
the energy and significantly to the specific heat of the 
normal state. We suggest that this is indeed the case for 
Tl 2201 (see below), for example. To resolve this ambi- 
guity, instead of the conjecture made by Anderson jlj, we 
propose to subtract the fluctuation effects and to use the 
remainder as an effective specific heat from which to ex- 
tract the c-axis contribution to the condensation energy. 
The rationale is that free energy can be decomposed into 
a singular and a non-singular part. The universal sin- 
gular part is more sensitive to collective long-wavelength 
fluctuations, while the non-singular part is dominated 
by short distance microscopic pairing correlations. This 
procedure is well suited to ILT, because, in this theory, 
the effective "mean field" condensation energy can be 
enhanced due to pair tunneling between layers . Note 
that there is no simple relation between T c and conden- 
sation energy, except in mean field theory. 

The fit to the specific heat of Tl 2201 to 2D Gaus- 
sian fluctuation plus non-singular terms 29|3(J is shown 
in Fig. 0. We have used C(T > T c ) = jT + g + /t 




20 40 60 80 100 120 140 T <K> 

FIG. 1. The electronic specific heat data of Tl 2201[15] 
fitted to a combination of singular and analytic terms (solid 
line); T c — 78.7 K. The straightline is 7T, and the dashed line 
is the analytic part of the specific heat below T c . 

and C(T < T c ) = c T(l + c x t + c 2 t 2 ) + g_/t, where 
t = |1 - T/T c \ (7 = 0.59, ff+ = 2.38,. g_ = 0.74, c = 
1.44, a = -2.79, c 2 = 2.07). The fit for 3D Gaussian 
fluctuations is considerably worse. This suggests that the 
in-plane contributions dominate. The data for Hg 1201 
are too imprecise to do the same, but clearly fluctua- 
tion effects are quite prominent and a proper subtraction 
should result in a larger prediction for A c . Optimally 
doped LSCO does not exhibit fluctuation effects that are 
as pronounced. For underdoped LSCO, we were unable 
to use the specific heat data |lj| as they do not seem to 
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fit any simple form. 

Condensation energies are obtained from an integra- 
tion of the measured specific heat. In Table Q we show 
both results with (AE su b) and without (AEo) subtract- 
ing fluctuation effects. By using these values for the right 
hand side of Eq. (||) we extract the corresponding values 
of A c as shown along with the experimental values [ p2[ . 
The precision of the present Tl 2201 specific heat data is 





LSCO (15%) 


Tl 2201 


Hg 1201 


AE 


~ 150 


~ 825 


~ 850 


A-E su b 




~ 25 




A c ,o 


~ 6.5 


~ 1.7 


~ 2 






~ 10 




Ac,exp 


5.5 H 


12 - 19 [2,5] 


6 - 8 [3, 12] 



TABLE I. Condensation energies (in mj/g-at) and pene- 
tration depths (in /im). Precise error estimates are unavail- 
able. 

not sufficient to make a precise subtraction of the fluc- 
tuation contribution. The uncertainty is considerable; a 
reasonable guess for the uncertainty in AE sub is 25 ± 15. 
Nonetheless, it is clear that using the subtracted value 
gives a much larger penetration depth. 

/. The Basov correlation: We manipulate the right 
hand side of Eq. (|^) to draw further conclusions. We 
perform a canonical transformation such that H c is elim- 
inated from the hamiltonian H . Thus. 



H = e~ s He s = H Test 



(5) 



where the antihermitian operator S is defined by H c + 
[H rest , S] = 0. The ground state 0) of the full hamilto- 
nian H can be determined perturbatively in S (or, equiv- 
alently t±) to show that the ground state expectation 
value of the H c is given by 



(6|# c |6) = -2^ 



KQIWI 

E„ — E 



(6) 



where E n and \n) are the eigenvalues and eigenfunctions 
of -ffrest- Of course, the same result could be obtained 
directly without making a canonical transformation. We 
have taken this route to hint that the canonical trans- 
formation, if carried out in infinitesimal steps, could po- 
tentially be a powerful method to obtain the effective 
low-energy hamiltonian p3| . 

For conserved parallel momentum, the expansion on 
the right hand side of Eq. (||) does not converge in a 
Fermi liquid theory because of vanishing energy denom- 
inators; therefore the expansion would not be valid. In 
a gapped state, the expansion can be legitimate because 
of the absence of vanishing energy denominators. In a 
non-Fermi liquid state, the matrix elements should van- 
ish for vanishing energy differences, and the the sum 



is skewed to high energies. Thus, the energy denom- 
inator can be approximated by W 0], and the sum 
can be collapsed using the completeness condition to 
(0|# c |6) w -(0|#f|0)/W. The effective hamiltonian 
— H 2 /W is identical to the hamiltonian of previous re- 
alizations of ILT ||,|6|. 

Thus (H c ) is of order t\/W. Then, from Eq. (§) for 
example, one can see that on dimensional grounds the 
c-axis conductivity is 



cr c (T) = a 



f e 2 dt\ \ 1 

\Awn 2 ) n(T)' 



(7) 



where a is a numerical constant weakly dependent on the 
band structure. The inelastic scattering rate is propor- 
tional to the unknown function f2(T). Combining the 
result of the previous paragraph with Eqs. (0,0), we find 



C 2 47T , , , , , 

^ = —a c {T)VL{T) [u s 



(8) 



where M Si „ is (0\(H c /t±) 2 \0) s ^ n . The average here is with 
respect to the ground state of Hrest) |0) not |0). The 
product o~ c (T)fl(T) is independent of T. 

For underdoped to optimally doped materials, the c- 
axis resistivity, p c (T), can often be fitted to 



p c (T) = b,T-P 



b' 2 T. 



(9) 



The logarithmic behavior |2(| obtains in the limit p — ■> 0. 
If we express Eq. (||) in terms of the temperature T* at 
which the resistivity takes its minimum value, we get 



c 

a! 



4tr7 c (T*)T* 



b 2 (p + l) 



P 



(10) 



where b 2 = b' 2 {de t ± /Ti AW). The expression in the 
curly brackets depends dominantly on b 2 , which describes 
the high temperature resisitivity. The low temperature 
behavior enters only through the exponent p. Thus, pro- 
vided the expression in curly brackets is a universal con- 
stant, a plot of In A c against \n[o- c (T*)T*] should be a 
universal straight line, independent of material, with a 
slope of —1/2. Basov et al. |[l| discovered a similar cor- 
relation by plotting lnA c against a c (T c ), shown as (I) of 
Fig. g, The correlation discussed here, shown as (II), is 
excellent. The data for underdoped LSCO, however, are 
affected by both the structural transition and the (1/8)- 
anomaly. In Fig. H, we have taken T* m T c for those 
optimally doped materials that show simply a flatten- 
ing of p(T) close to T c . Thus, we see that b 2 is indeed 
inversely proportional to (u s —u n ), which, in ILT, is pro- 
portional to the T = superfluid density, n s (0). This 



can be tested further in future experiments |33|. 

g. Conclusion: ILT accounts for a number of exper- 
imental behaviors, in particular the Basov correlation, 
and it provides a recipe for determining the c-axis pen- 
etration depth. In Tl 2201 and Hg 1201 there must be 
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I) 



LSCO 1 



Hg 1 
T 



Hg2 
T 



a ( T c ) (ohm" 1 cm" 1 ), T*o( T*) (Kohm" 1 cm" 1 ) 

FIG. 2. The Basov plot: lnA c is plotted against \na c (T c ), 
as in the original Basov plot (I), and against ln[T* a c (T*)] as 
discussed here (II) . The legends in group (II) are the same as 
those in group (I). Til: Ref. [2]; T12: Ref. [5]; Hgl: Ref. [3], 
Hg2: Ref. [13]; LSCO 1 (12%), LSCO 2 (15%): Refs. [15, 28, 
32]. 



strong superconducting correlations in the normal state. 
The source of these must be both the fluctuation effects 
not contained in the mean field treatment of ILT as well 
as substantial in-plane pairing correlations. Although 
ILT is not the main driving mechanism for Tl 2201, it 
may be for LSCO and in any case ILT remains an im- 
portant mechanism which can enhance T c in both single 
layer and multilayer materials [|7| ]To| J25[ | . 
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